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■ Abstract. I briefly review the theoretical models of radiatively efficient, 

geometrically thin and optically thick accretion disk spectra that currently exist 
. for AGN. I then discuss three recent observational developments that have real 

potential to teach us about the physics of these flows. Finally, I present results on 
the most recent, thermodynamically consistent simulations of magnetorotational 
, turbulence and discuss what these simulations are suggesting about the vertical 
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structure of accretion disks 



> 

OO . 1- Current Standard Accretion Disk Models of AGN 

m 

Accretion disks around supermassive black holes continue to be the generally 
accepted paradigm for the central engine of AGN. I focus in this review on "stan- 
dard disks" : those that are geometrically thin and optically thick. These models 
are, by construction, radiatively efficient. For radiatively inefficient disks, see the 
papers by Mineshige and by Yuan in these proceedings. Currently all theoretical 
models of the photon emission of accretion disks (radiatively efficient or other- 
wis e) are based on the alpha prescription for the anomalous stress introduced 
by IShakura Sunvaevl (jl973l ) . A large numbe r of such mod e ls for standard 
disks h ave been constructed o yer the years, e.g. Mal kanl (11983T ) ; ISun fc Malkail 



(|1989T ) ; lLaor & Netzerl (|1989T ) ; lLaor. Netzer Piranl (|l990l ). to name just a few 



The most sophisticated of thes e mod els h ave been developed in a ser ies of 
rS papers bv lHubenv fc Hubenvl (|l997t Il998l ) and lHubenv et al.1 fcoOOl . l200ll ). and 



are available to the community through the IDL package AGNSPEC, written 
by Ivan Hubeny. The models are constructed by first dividing the accretion 
disk into a set of discrete annuli, and then computing a one-dimensional vertical 
interior and atmosphere structure based on certain assumptions in each annulus. 
Using the locally calculated spectrum, photon geodesies from the disk to infinity 
are computed to then produce a predicted model spectrum for an observer at 
some viewing angle. Some aspects of the physics of accretion disks are treated 
very well by the AGNSPEC models. They ful l y account for the relativistic 
radial structure of the disk (iNovikov Thornd Il973l : iRiffert & Heroldl Il995l ) 



and relativ istic Doppl er shifts, gravitational redshifts and light bending in a Kerr 
spacetime (Ago! 19971 ). The models include a detailed non-LTE treatment of all 



abundant elements and ions in the disk, and include continuum opacities due 
to bound-free and free-free transitions of these ions, as well as Comptonization. 
(No lines are included in the models, however.) 

There are many ad hoc assumptions that continue to plague these models, 
however. First there are assumptions that relate to the radial structure of the 
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disk. The disk is assumed to be stationary, with a zero torque inner boundary 
condition at the innermost stable circular orbit. Recent global simulations of 
MRI turbulent flows around black holes suggest instead that substantial torques 
on the disk can be exerted by material in the plungin g region through magnetic 
field lines that connect this mate rial to the disk (e.g. Krolik. Hawlev. Hirosel 
20051 ). The models also adopt the lShakura &: Sunvaevl ( 19731 ) alpha prescription 



for the anomalous stress, with the stress being proportional to the total ther- 
mal (gas plus radiation) pressure and the alpha parameter assumed constant 
with radius. This is done even though such models, if taken literally, are ther- 
mally and "viscously" uns table when radiation pressure dominates gas pressure 
( Shakura Sunvaev| [i976). Reprocessing of radiation that leaves the disk at one 
point and returns to a different point due, e.g. to light bending, is neglected. 
This can be an im portant effect particularly for rapidly spinning black holes 
( Cunningham! [19761 ) . 

There are perhaps even more serious assumptions regarding the vertical 
structure at any radius. This vertical structure is first assumed to depend only 
on height, i.e. it is entirely one dimensional. As I discuss below, this is prob- 
ably not accurate as substantial inhomogeneities probably exist. The structure 
is also assumed to be symmetric about the equatorial plane of the black hole: 
no warps or bending waves are considered, nor are time-dependent asymm etries 
induced by the turbulence accounted for ( Hirose. Krolik. Stond [20061 ) . The 
disk is assumed to be in vertical hydrostatic equilbrium with gas and radiation 
pressure being the only means of support against the vertical tidal field of the 
black hole. Magnetic forces are neglected. Just as a model of a stellar interior 
requires specification of the sources of nuclear energy, a model of the vertical 
structure of the disk requires specification of the vertical profile of turbulent 
dissipation. The AGNSPEC models assume that the dissipation per unit vol- 
ume tracks the density, i.e. the dissipation per unit mass is constant. Heat 
is also assumed to be transported vertically through radiative diffusion. These 
assumptions are not consistent in a radiation pressure dominated annulus: con- 
stant dissipation per unit mass plus radiative diffusion immediately leads to 
a density which is constant with height if the opacity is dominated by elec- 
tron scattering (a good approximatio n). This in turn is convectively unstable 
(jBisnovatvi-Kogan &: B linnikovl Il977l ) . 

It is probably for many of these reasons that the model s do not do a very 
good job in reproducing some aspects of the observations (e.g. iKoratkar Blaea 
1999). For example, they do not produce X-rays. They prod uce a break at th e 
Lyman limit of hydrogen but not quite like what is observed ( Blaes et al.ll200il ). 
nor do the expected trend s in far ultraviolet continuum slopes agree with obser- 
vation ( Shang et al.| [2005). The models cannot explain th e observed near simul- 
taneo us optical/ultraviolet variability ( Alloin et al.lll985l : ICourvoisier Clavel 
199ll ) because they do not include reprocessing. 

In spite of the fact that the physics of these models is so uncertain, I should 
emphasize that exactly the sam e type of modeling app lied to stellar mass black 
hole accretion disks (BHSPEC. iDavis &: Hubenv 20061 ) results in very good fits 
to the spectra of black hole X-ra y binaries in the high-soft/thermal dominant 
state (jDavis. Done. &: Blaesl 120061 ) . For LMC X-3, for example, where we know 
the distance and the black hole mass, we can fit the accretion disk spectrum very 
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well over a range of luminosities with fixed spin and inclination. The models only 
fail at the very highest luminosities (Eddington ratio greater than of order 0.7), 
where the assumption of a geometrically thin disk is almost certainly wrong. 
Either the ad hoc physical assumptions that go into the disk models are not too 
wrong, or the spectra are robust to changing these assumptions. 



2. Some Recent Observational Developments that Provide New Con- 
straints on Disk Physics 

I would now like to mention three recent observational developments that are 
exciting because they provide us with new empirical constraints on what is going 
on in AGN accretion disks. 

2.1. Spectropolarimetry 

For good reason, the broad and narrow emission line regions in AGN have 
attracted enormous attention from the community for years. However, to a 
parochial theorist like myself who is primarily interested in the accretion flow 
very near the black hole, emission from these regions merely contaminates our 
view of the innermost parts of the central engine. In particular, the numerous 
broad and narrow emission lines, the small blue bump, and the infrared thermal 
emission from dust make an accurate determination of the shape of the big blue 
bump (which is supposed to originate in the accretion disk) very difficult. 

Fortunately, all this contaminating emission appears to be unpolarized in a 
subset of quasars, whereas the underlying big blue bump is polarized. The 
reasons for this are unclear, but it appears that there is an electron scat- 
tering medium enveloping the central accretion flow but lying entirely within 
the broad line region in these quasars. If one measures the polarized flux 
in these quasars, all the contaminating emission is removed and the shape 
of the underlying continuum is revealed. In particular, an absorption feature 
which is almost certainly a Balmer edge lies underneath the small blue bump 
(jKishimoto. Antonucci. Blaed 120031 ; iKishimoto et al.ll20()3 ). This is the first 



direct evidence that the big blue bump originates from thermal emission from an 
optically thick medium! This technique has also been used to measure the long 
wavelength, near-in frared shape of the big blue bump un derneath the thermal 
emission from dust ( Kishimoto. Antonucci. Sz Blaes|[2005l ). Kishimoto discusses 



these results in much more detail in his paper in these proceedings. 
2.2. Microlensing 

While it might be possible to directly image the very central part of the accretion 
flow around the bla ck holes in the Galactic center and M87 in the not so distan t 
future with VLBI (|Falcke. Melia. Agoll l2000l : iBiretta. Junor. Livid [2002). 



the angular sizes of the radiatively efficient black hole accretion disks that we 
know about in both AGN and X-ray binaries are hopelessly small for imaging to 
really be practical. Fortunately, nature has provided us with a cosmic magnifying 
glass in the form of gravitational lensing. Many examples of multiply imaged 
quasars are now known and are fit reasonably well by models of the foreground 
lensing galaxy (the macrolens). Depending on the source size at a particular 
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observed wavelength, each of the quasar images can in turn be microlensed by 
stars that reside in the foreground macrolens galaxy. If the source (the quasar) 
has a size scale much less than the Einstein radius of a typical stellar microlens, 
then the flux ratios of the different quasar images will typically differ markedly 
from the predictions of the macrolens model. They will also vary with time in a 
manner that has nothing to do with the intrinsic variability of the quasar. On the 
other hand, if the source is extended compared to the microlens Einstein radius, 
there will be little differential magnification and the flux ratios should agree with 
the macrolens model. Microlensing therefore probes the size of quasar emission 
regions on scales comparable to the Einstein radius of a stellar microlens. For 
cosmologically distant lenses and sources, this corresponds to angular resolutions 
of several microarcseconds, or linear scales in the quasar ~ 5 x 10 16 cm, for a 

sola r mass microlens. 

Rauch fc Blandfordl (|l99ll ) were the first to use this technique. For the Ein- 



stein Cross, the qu adruply lensed quasar where microlensing was first discovered 
(jlrwin et al.lll989l ). they found that the optical (really near ultraviolet) emission 
regions were smaller than expected for thermal emission from an accretion disk. 
More modern work in this system and other gravitational lenses has generally 
reached the opposite conclusion: the near ultraviolet emission regions are larger 
than predicted by simple disk models. In particular, iPoolev et al.l ( 20061 ) have 
found that near ultraviolet half-light radii predicted by multitemperature black- 
body disks are a factor 3 to 30 times too small than what is necessary to explain 
the observed microlensing flux ratios in ten quasars, including the Einstein Cross. 
Using a different analysis of the microlensing stat istics and (hopefully b etter) 
black hole mass es based on emission line widths, Kochanek et al. (2006) and 



Dai et al.1 (120061 ) also find that the near ultraviolet emission regions are larger 



than would be expected for the size of the annulus with the corresponding effec- 
tive temperature. However, they also find that the size of the emission region 
has roughly the right scaling with black hole mass betw een different sou rces. 

It is not yet clear what causes the discrepancies ( Dai et ail 120061 ) . Con- 
tamination of the observed fluxes by the small blue bump is clearly a worry, 
as this originates from larger scales than the disk. Reprocessing of inner disk 
radiation will also increase the size of the near ultraviolet emitting annuli. Disk 
atmosphere models such as AGNSPEC show significant deviations from local 
blackbody emission: the theoretical spectrum is harder, which would make the 
size of the emission region at a given wavelength larger than expected from a 
local blackbody disk. Whatever the cause, this observation technique provides a 
unique probe of the size scales of the emission regions, and I strongly encourage 
observers to extend it across multiple wavelengths. 



2.3. SDSS Spectra with Black Hole Masses 

As discussed by Kaspi and Peterson in these proceedings, reverberation mapping 
leveraged by the broad line region radius/continuuum luminosity correlations 
has given us a simple method to measure approximate black hole masses, and 
this can be done for huge numbers of quasars in the Sloan Digital Sky Survey 
(SDSS) sample. As a result, one can now look at continuum colors as a function 
of luminosity L and black hole mass M and see whether the variations match 
the trends predicted by accretion disk theory. 
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Bonning et al. ( 20061 ) have taken a first stab at this by looking at col- 



ors at three pairs of rest frame wavelengths (5100A/4000A, 4000A/2200A, and 
2200A/1350A) versus a measure of the expected maximum disk effective tem- 
perature T max oc L^M -1 / 2 . They found that the 5100A/4000A colors pre- 
dicted by the AGNSPEC models agreed with the SDSS data at all T max . The 
4000A/2200A colors agree at low values of T max , but not at high, and the 
2200A/1350A colors agree only at the lowest values of T max . The sense of the 
disagreement in the discrepant regions is that the data is redder than the AGN- 
SPEC predictions. Internal reddening by dust in the quasars is therefore a 
worry, but they argue that this cannot explain the observed discrepancies. In- 
terestingly, they find that the value of T max which marks the point where the 
observed 4000A/2200A colors depart from the model predictions happens to cor- 
respond to an Eddington ratio of 0.3, beyond which the standard geometrically 
thin disk assumption may break down. "Slim" disk physics may therefore be 
playing a role (see the paper by Mineshige in these proceedings) . 

As iBonning et al. ( 20061 ) acknowledge in their paper, the combination of 
luminosity and mass that appears in T max actually depends empirically only on 
the measured width of the Mg II line used in their analysis, and not on the 
measured luminosity. It would be nice to extend their work to a truly two- 
dimensional analysis with respect to luminosity and black hole mass. 



3. Some Recent Theoretical Developments: a Begone! 



The IShakura &: Sunyaevi ( 1973! ) alpha prescription that has plagued accretion 



disk models for decades is slowly being banished to oblivion. Thanks to th e dis- 
covery of the magnetorotational instability (MRI. iBalbus Hawlev 1991 ) and 



the ever increasing power of numerical simulations, the theoretical community is 
starting to rebuild accretion disk theory from scratch. One of the most exciting 
developments is that radiation magnetohydrodynamic simulations are telling us 
how the vertical structure might really work. Such simulations are possible in 
a local vertical slice through a disk, including the vertical gravity and the MRI 
turbulence, capturing the dissipation of the turbulence as heat, and allowing 
radiation to diffuse out of the box so that cooling can occur. Both the dy- 
namics and the thermodynamics of the turbulence are therefore captured in the 
simulations. 

Three such simulations have been completed so far. The first was by I Turner 



( 20041 ) . who modelled a radiation pressure dominated annulus and found that 



it reached a thermal equilibrium, albeit with large fluctuations. There was no 
evidence of a thermal instability, but there was significant vertical expansion of 
the disk during a heating phase which caused approximately half the mass of the 
annulus to leave the grid. The dissipation per unit volume inside the annulus 
did not track the density, and the average vertical structure was convectively 
stable. Time step constraints on the radiation diffusion solver required the use 
of a rather high density floor within the simulation which meant that every zone 
was optically thick: there was no photosphere within the simulation domain. 
The results of this pioneering simulation are tantalizing, but they are not yet 
definitive. 
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A number of improvements have since been made to the code used by I Turner 

(2004) so that it conserves energy much better and also runs more efficiently and 
faster, allowing the density floor to be reduced to a level where the photosphere 
can now be included within the simulation domain. A truly radiation pressure 
dominated annulus has not bee n successfully run yet with t he new code, however, 
possibly because the standard IShakura Sunvaevi ( 19731 ) vertical structure as- 
sumptions are so far from reality that using them as an initial condition leads to 
very larg e transients. A simulation of a g as pressure dominated annulus was pub- 
lished bv lHiroseT Krolik. & Stond" (|200a ). A new simulation with co mparable gas 
and radiation pressures a t the midplane has now been completed ( Blaes et al.l 
2006bl : iKrolik et al. 2006). It is probably the best simulation we have so far of 



the conditions that might exist in black hole accretion disks, although it is still 
not radiation dominated. 

Figure 1 sh ows the time variat ion of the energy content in the box for this 
new simulation (|Krolik et al.ll2006h . The average cooling time in the annulus 
is only ~ 7 orbits, and clearly the energy content of the box goes through 
very large fluctuations on time scales much longer than this. The nominal 
Shakura Sunvaev ( 19761 ) thermal instability criterion that total radiation pres- 



sure exceed three fifths of the total thermal pressure is violated around 90 orbits 
and also (slightly) around 300 orbits. Yet there is no evidence of a thermal 
runaway. The reason for this appears to be that the volume-integrated dissipa- 
tion rate is approximately proportional to the square root of the total energy 
content during the hottest epochs, while the cooling rate is proportional to the 
total energy content. These scalings therefore lead to thermal stability. It is far 
from clear that this will extend to regimes of much higher radiation pressure, 
however. 

That describes the thermal equilibrium, but the hydrostatic equilibrium is 
also interesting. Away from the midplane, vertical magnetic forces dominate gas 
and radiation pressure gradients in supporting the annul us against gravity. This 
was al so seen in the gas pressure dominated simulation bv lHiroseT Krolik. & Stond 
(2006). This has immediate spectral implications as it produces a larger density 
scale height compared to what would exist if only gas and radiation pressure 
gradients were taken into account. This results in a lower density at the effec- 
tive photosphere, producing increased ionization and red uced absorption opa city. 
This generically results in a harder emergent spectrum ( Blaes et al.ll2006al ). 

However, one-dimensional ste llar atmosphere modeling does not do the sim- 
ulation justice ( Blaes et al. 2006b). As shown in Figure 2, strong horizontal den- 
sity fluctuations (up to a factor 100) exist as deep as the effective photosphere of 
the horizontally averaged structure and even deeper. Some of the high density 
regions are reminiscent of the shock train st ructure predicte d in the nonlinear 
development of the photon bubble instability ( Begelman1l200ll ) . but the expected 
correlations between radiation flux and fluid velocity do not appear to match 
this instability. On the other hand, the transverse Parker instability is clearly 
present: high (low) density regions at depth are strongly correlated with upward 
(downward) magnetic tension forces. 

At least to the extent that thermalization of photons with matter occurs 
through true absorption and emission processes (as opposed, say, to Compton 
scattering) the large density fluctuations will tend to enhance the thermaliza- 
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Figure 1. Time variation of the total energy content for the new simula- 
tion with comparable gas and radiation pressures. From top to bottom, the 
curves represent the vertically integrated radiation energy, gas internal energy, 
magnetic energy, and fluid kinetic energy, respectively. 



tion of photons, leading to a more blackbody-like, softer emergent spectrum 
(jDavis et al.ll2004l : lBegelmadl2006h . 

The complex photosphere will also presumably reduce the degree of polar- 
ization of the emergent radiation, as the a tmosphere no longer has the plane- 
parallel symmetry used in simple models ( Chandrasekhai I960] ). This has of 
course been a long-standing problem with disk models: electron scattering dom- 
inated slabs gen erally produce too high a polarization which is also at the wrong 
position angle ( Stockman. Angel. &: Milev 19791 ; Antonucci 1988). The inho- 
mogeneities that we see in the simulation resemble the "rough surface" that 
Coleman fc Shields! (119901) proposed might reduce the polarization. As first sug- 
gested by Gnedin Sz Silant'evl (jl978l ). Faraday depolarization might also be sig- 
nificant. The Faraday rotation angle at the peak of the locally emitted thermal 
spectrum is roug hly 0.8t Tj R 1/2 , where tt ~ 1 is the Thomson depth of the pho- 
tosphere and R is the ratio of magnetic energy density to radiation pressure at 
the photosphere. For the conditions observed in the simulation at times of high 
energy content, this angle turns out to be approximately 2 radians. 

To summarize, the simulation suggests that the vertical structure of a disk 
annulus with comparable gas and radiation pressures is as follows. MRI turbu- 
lence churns away at the midplane where the vast majority of the accretion power 
is being dissipated as heat. Away from the midplane, but still well inside the 
photosphere, magnetic energy density dominates both gas and radiation pres- 



8 



Blaes 



i°g, n [p/(g cm 3 )] 



-6.8 



-6.1 



4.6 



-3.9 



-3.2 



-2.4 



4.5 



4.0 



N 



3.5 - 



3.0 



2.5 




Figure 2. Density (gray scale) and projected magnetic field (arrows) in a 
constant radial slice through the subphotospheric upper layers of the sim- 
ulation at 90 orbits. The upper and lower horizontal dashed lines show the 
Rosseland mean and effective (i.e. thermalization) photospheres, respectively, 
of the horizontally averaged structure. 



sures, and magnetic forces provide the hydrostatic support of the disk against 
gravity. These magnetically supported layers are Parker unstable, resulting in 
a very complex, inhomogeneous photosphere higher up. Outward vertical heat 
transport is completely dominated by radiative diffusion, and not by advection, 
convection, or Poynting flux. Because most of the dissipation is occurring at the 
midplane, no energetically significant corona is produced above the photosphere. 

The absence of a corona was al so seen in the gas-pressure dominated simula- 
tion ( Hirose. Krolik. &: Stone! 20061 ) . and is a potential problem for models which 
suggest that much of the accretion power can be carried locally out of the dis k by 
magnetic fields (e.g. Haardt &: Maraschi 199ll ; ISvensson & Zdziarskilll994T ). It 
might be that the simulation box is too small to allow enough magnetic field to 
get out and dissipate to produce a corona. (It turns out that the wavelength of 
the fastest growing Parker mode barely fits inside the box.) However, this may 
not change the fact that most of the dissipation occurs in the MRI turbulent 
regions near the midplane. It is important to bear in mind, though, that this 
dissipation is mostly numerical in nature, being caused by losses of magnetic 
and kinetic energy at the grid scale. Still this may be okay if MRI turbulence 
always cascades energy down from large to small length scales. 

If the simulation results are correct, then where do the X-ray emitting re- 
gions in AGN come from? It may be that things change when radiation pressure 
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dominates gas pressure at the midplane. Mag netic buoyancy may be enh anced 
in radiation pressure dominated environments ( Sakimoto Sz Coroniti 1989), and 
perhaps this would alter the vertical dissipation profile to put more power out 
in a corona. Ho wever, this was not observed in the preliminary simulation 
by iTurnerl (|2004ft . It could also be that the X -ray source has nothing at all 
to do with the local disk, but is instead material originating from the inner- 
most disk and outflowing in a jet or backflowing over the rest of the disk, 
as is seen in global, n onradiative, general relativistic MHD simulations (e.g. 
Hawlev k Balbusl I2002T ) . Regions of high electric curre nt density are also seen 



near and inside the plunging region in such simulations ( Hirose et al.ll2003 ). and 
the dissipation associated with such high currents might very well produce an 
X-ray emitting coronal plasma. 

I hope I have given the reader a taste of the progress we have made in under- 
standing the fundamental physical processes inside accretion disks. Radiation 
pressure dominated simulations are now in the works! 
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